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ABSTRACT 

Context. It has been demonstrated that active galactic nuclei are powered by gas accretion onto supermassive black holes located at 
their centres. The paradigm that the nuclei of inactive galaxies are also occupied by black holes was predicted long ago by theory. 
In the last decade, this conjecture was confirmed by the discovery of giant-amplitude, non-recurrent X-ray flares from such inactive 
galaxies and explained in terms of outburst radiation from stars tidally disrupted by a dormant supermassive black hole at the nuclei 
of those galaxies. 

Aims. Due to the scarcity of detected tidal disruption events, the confirmation and follow-up of each new candidate is needed to 
strengthen the theory through observational data, as well as to shed new light on the characteristics of this type of events. 
Methods. Two tidal disruption candidates have been detected with XMM-Newton during slew observations. Optical and X-ray follow- 
up, post-outburst observations were performed on these highly variable objects in order to further study their classification and 
temporal evolution. 

Results. We show that the detected low-state X-ray emission for these two candidates has properties such that it must still be related 
to the flare. The X-ray luminosity of the objects decreases according to theoretical predictions for tidal disruption events. At present, 
optical spectra of the sources do not present any evident signature of the disruption event. In addition, the tidal disruption rate as 
derived from the XMM-Newton slew survey has been computed and agrees with previous studies. 
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1 . Introduction much debate concerning LINERs and their powering mechanism 

. 5^ caused either by an active nucleus or starburst. Given their activ- 

k> ; According to our current picture of galaxy demography, massive j^ jgyel and spectral properties they could represent the connec- 

*\1 . dark objects could reside at the centres of many galaxy bulges tjo^ between active and inactive galaxies ( Gonzalez-Martin et al| 

c3 ■ JKormendv & Richstonelll995H Magorria n et alJll998 ^ While it f2(Ml . 

■ - - ' has been demonstrated that active nuclei powered by gas accre- ^h^ paradigm that the cores of inactive galaxies are occupied 

tion onto supermassive black holes (SMBHs) populate the cores ^y quiescent SMBHs was predicted long ago. The discovery of 

of a number of galaxies, confirmation of their dormant state lurk- dynamical evidence of a mass concentration with properties of 

ing in non-active galaxies is difficult to obtain. j^is type of object at the centre of our own galaxy and after- 

Chronologically, quasars were more abundant in the early wards in a handful of o ther galaxies like M31 and N GC 4258 

Universe at z~2 than at present so dead quasar engines are ex- favours this conjecture JKormendv & RichstonellT995 '). An un- 

pected to be enclosed in the nuclei of otherwise inactive galax- avoidable consequence of the existence of remnant SMBHs at 

ies. Evidence connecting the genesis of SMBHs and the pos- the nuclei of optically inactive galaxies is the detection of flare 

sible hierarchical quasar to galaxy evolution remains elusive, radiation pro duced when a star is tida l ly disrupted and accreted 

At the lower end of the activity scale and representative of the by SMBH s (iLidskii & Ozernoilll979l: ICarter &Luminetlll982l: 

most common activity level in the local Universe r eside low- |Reesll988l) . These occurrences, so-called tidal disruption events. 

ioniza tion nuclear emission-line regions (LINERs) (|Heckman| would generate outbursts of EUV/X-ray radiation decaying on a 

|1980|). Their optical spe ctra are similar to those of giant ellipti- time scale of several months to years. Although active galactic 

cal galaxies dPenicoloet al. 2005 ), dominated by low-ionization nuclei (AGN) can also be susceptible to these phenomena, the 

emission lines, that are also similar to fl ie iiarrow lines in most unambiguous cases for a stellar disruption come from host 

Seyfert galaxies but with different line ratios (|Ho|2008|). There is galaxies with no - or only faint - permanent activity. 

It is expected that tidal flares emit a luminosity peaking in 

Send offprint requests to: P. Esquej the EUV/X-ray regime. When assuming a black body radiat- 
e-mail: pili@mpe.mpg.de ing at the Eddington luminosity, the temperature of the emit- 
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ting region calculated at the last stable orbit (3rs, where rs is the 
Schwarzschild radius) will be ~60 eV for a 10^ Mq black hole. 
At low redshifts, the optical and UV bands will be sensitive to the 
Rayleigh-Jeans tail of the black body radiation, but these events 
will be best detected in X-rays, with giant amplitudes of vari- 
abih ty that will reach factors of hundreds up to several thousand 
(e.g.lKomossa et al.l l2004HHalpern et ai]|2004l) . 

The comparison o f sources from th e ROSAT PSPC All- 
Sky Survey (RASS) (IVoees et all fT999l) with ROSAT PSPC 
pointed observations allowed the detection of five tidal dis- 
ruption candidates. These sources showed extreme X-ray vari- 
ability and properties in agreement with those expected for 
tidal disruption events (see Ko mossa.2002^ for a review): ex- 
treme X-ray softness in outburst. X-ray peak luminosity up 
to ~ 10'*^ erg s~' that declines at a rate r^^^ on a timescale 
from months to year s ([Bade et al. 1996; K omossa & B ade 1999; 
iKomossa & Greineii [l999l; iGrupeetalJ [1 299|). A few candi- 
dates have emerged in the Uy/optical (e.g.lRenzini et aill 19951; 



IStern et al.ll2004FlGezari et"aLll2ob8HKomossa et al.ll2008l) . three 
of which were subsequently detected in soft X-rays between 1 
and 3 years after the high-state detection. 

These occurrences are very rare in the sky. Theoretical mod- 
els postulate that a tidal disru ption event sh ould occur once ev- 
ery 10^ - 10^ yr per galaxy dWang & Me rrit 2004) depending 
on the stellar density in the nuclear cusp and the SMBH mass. 
The derived tidal disruption rate from the RASS, although in- 
cludi ng large uncertainty, lies in agreement wi th these predic- 
tions (ISembav & Westlll993l:lDonlev et al.ll2002 ). Although pre- 
viously detected tidal disruption events support theoretical pre- 
dictions, due to the paucity of candidates the detection and con- 
firmation of each new case will shed new light on the phenomena 
here debated. 

The identification of tidal disruption events discussed in 
the present paper presupposes the existence of two large area 
X-ray sky surveys performed at different epochs. Correlating 
the XMM-Newton Slev/ Survey Source Catalogue (XMMSLl) 
JSaxton et al.l2008l) with the RASS we foun d two galaxies whic h 
were candidates for tidal disruption events ("Esquej et al.ll2007l) . 
Here we present follow-up multiwavelength observations per- 
formed to investigate the classification of these sources and 
study their temporal evolution. 

This paper is structured as follows: in section (|2]i we intro- 
duce the basics of tidal disruption phenomena, then we define 
the criteria followed for selecting our sample and present the X- 
ray and optical properties of the sources. Detailed analysis of 
follow-up observations is described in section (|4]l. X-ray light 
curves and estimations of the black hole masses are calculated 
in the subsequent sections. The tidal disruption rate as derived 
from the XMM-Newton slew survey is derived in section (|7J and 
our conclusions are finally summarised. 

Throughout this paper we will assume a ACDM cosmology 
with (Qm, ^a) = (0.3, 0.7) and Hq = 70 km s"' Mpc '. All X- 
ray fluxes and luminosities here quoted have been corrected for 
Galactic absorption. 



2. Theory of tidal disruption events 

A star of mass m* and radius /?* orbiting around a SMBH of 
mass Mbh will be disrupted when it approaches the black hole 
tidal radius 



as this is the distance where the surface gravity of the star equals 
the tidal force from the hole (yU is a dimensionless coefficient of 
order unity). The stellar disruption process is predicted to happen 
up to black hole masses of 10*^ Mq. For objects above this limit, 
the tidal radius is inside the Schwarzchild radius and, therefore, 
a solar-type star will fall into the black hole without disruption. 
Atmosphere stripping of giants is still possible for larger black 
hole masses. 

During its infall towards the central black hole the star ex- 
periences a gravitational torque which spins it up and creates a 
large velocity dispersion within the disturbed material. Once the 
star is disrupted, half of its debris will have enough energy to 
exceed the escape velocity and will be consequently ejected on 
hyperbolic orbits. The remaini ng half will be bound, ret urning 
to pericentre and circularizing f UlmeiJ[l999l; lReeslll988b . This 
returning material that initially follows highly elliptical orbits, 
will form an orbiting torus around the black hole due to dissipa- 
tion and internal shocks. The observable flare of radiation will 
begin at a time coincident with the accretion of the most tightly 
bound material as it first returns to pericentre. After one post- 
disruption orbit, a fraction of the remaining bound debris will 
be accreted by the hole at about the Eddington rate. lAval et al] 
(2000) showed that as little as 10% of the total mass of the star 
would be accreted, unlike previous estimates that claimed that 
all the bound material (50% of the star) would be accreted by 
the black hole. During this stage the soft X-ray luminosity of 
the galaxy is increased by up to several orders of magnitude and 
then declines at a rate t^^^ ^ on a timesca le from months to years 
(lEvans & Kochanekll989tlUkneiill999l) . 

The tidal disruption rate has been predicted by several au- 
thors based on the standard model o f collisional los s cone 
repopulation in a s p herical nucleus (iMagorrian et al.l 119981 ; 
ISver&Ulmel 119991: IWang & Mert3 |2004|). In general, this 
quantity depends on the mass of the black hole, the stellar den- 
sity, velocities in the galactic nucleus, and the rate at which radial 
loss cone orbits are depleted/replenished. Wang & Merrit (2004) 
derived a tidal rate of 10 "* - 10"^ yr ' Mpc"^ for a 10^ Mq black 
hole. 



3. Candidates selection 



^T = ^iR* 



Mr 



1/3 



(1) 



The state-of-the-art satellite XMM-Newton (iJansen et al.ll2001l) 
features high spectral resolution and excellent sensitivity due to 
the great collecting area of its mirrors coupled with the high 
quantum efficiency of the EPIC detectors. X-ray objects detected 
with the EPlC-pn camera on board XMM-Newton while the 
satellite is maneuvering between observation targets comprise 
the XMMSLl source catalogue. 

Our sample selection criteria aimed at finding tidal disrup- 
tion candid ates through comp arison of the XMMSLl with the 
RASS (see lEsquei etal.ll2007l for details). For the two result- 
ing objects, NGC 3599 and SDSS J132341.97-H482701.3, no 
prior AGN-like activity had been recorded from optical obser- 
vations according to previous published measurements (but see 
section [ni l. Both sources showed high variability and a soft X- 
ray spectrum as regards the slew observations. All these prop- 
erties favoured the interpretation of the flaring events as tidal 
disruption candidates. 

3.1. NGC 3599 

The slew source XMMSLl Jl 11527.3+180638 was detected 
in 2003-11-22 with an EPlC-pn 0.2-12 keV count rate of 
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Fig. 1. Detailed fits to the optical spectrum of NGC 3599. Flux calibrated data from ' Caldwell et aTl (l2003h are presented in black. 
The best fitting stellar population as derived from the EL-ICA approach is shown in red. The subtracted spectrum in blue illustrate 
the galactic nuclear region. Line identifications are given in the figure. 



5.5 cts s '. Through catalogue cross-correlation it was found 
that its counterpart is the SAO galaxy NGC 3599 at redshift 
0.0028 (iDenicolo et all l2005h . Within the early-type galaxies, 
NGC 3599 is a power-law galaxy (its brightness profile has a 
power-law shape) with a compact nuclear component not asso- 
ciated with non-the rmal acti vity (like those found at the centres 
of M31 and M33) dLauer et al. 1995). Although no public op- 
tical spectroscopic observations of this source were available, a 
detailed literature search revealed that NGC 3599 was already 
briefly included in other samples with different aims. Flux cali- 
br ated optical dat a were kindly provided to us from the sample 
of lCaldwell et al.l ( 12003) and then used to characterise the source 
in the pre-outburst stage. 

The long-slit optical spectrum of NGC 3599 was taken in 
2000-02-05 at the F. L. Whipple Observatory with the FAST 
spectrograph located in the 1.5m Tillinghast telescope (see 
ICaldwell et al.l I200 J, for observational details and data reduc- 
tion procedure). Prominent absorption features in combination 
with emission lines appear in the integrated spectrum. In order to 
further study the emission-line spectrum we first subtracted the 
starlight component. This step in the data analysis reduction is 
crucial to derive an accurate nuclear source classification based 
on emission lines ratios. Two independent approaches to the stel- 
lar population reconstruction have been implemented both using 
high resolution synthetic spectra from the evolutionary synthe- 
sis models of Bruzual & Chariot (2003) and assuming a Salpeter 
initial mass function (Salpeter 1955). These techniques comprise 
the Ensemble Learning Independent Component Analysis (EL- 
ICA) and the CONFIT codes. Discussion about technical details 
on both resources is beyond the scope of the current paper and 
presented elsewhere in the literature (see references below). 

The EL-ICA code presupposes that the stellar contribution 
is a linear combination of six non-negative independent com- 
ponents (ICs) that arise fro m the compaction of the 13 2 6 simple 
stella r populations (SSPs) of lBruzual & Chariot! ( l2003h dLu et alj 
l2006l and references therein). Through minimization of the x^^^ 



parameter, the optimum fit to our data is obtained with a single 
stellar population of 2.1 Gy and solar metallicity. Fig. [1] shows 
the data (black), stellar modelling (red) and the subtracted spec- 
trum that accounts for the galactic nuclear emission (blue). 

Conversely, the CONFIT approach consists of modelling 
the continuum spectral shape using two stellar components 
allowing them to have different redde ning and apply i ng y^ 
statistics to derive the opt imum fit dRobinson et al.l 120001 : 
Rodriguez Za urin et al.ll2007h . The continuum flux is measured 
in a number of wavelength bins (54 in this case) distributed over 
the spectral regions free of emission lines and atmospheric ab- 
sorption features. By computing the contribution of the different 
stellar components to reproduce the complete spectrum, the best 
fit resulted from a single stellar population of 2 Gy, which is 
fairly consistent with the outcome of the EL-ICA approach. . 

Spectral emission lines have been fitted using a single- 
component Gaussian. NGC 3599 exhibits narrow emis- 
sion lines of about 100 km s"' width. Intensity ratios of 
[O III] 5007 A/H/3=3.5(+ 1.0), [N II] 6584 A/Ha=l.l(+ 0.2) 
were measured from the analysis of the decomposed spectrum as 
returned by the EL-ICA code. As a comparison, those calculated 
from the CONFIT procedure provide [O III]5007 A/H/3^3.0(± 
0.2) and [N II]6584 A/Ha^l.li+ 0.3). Both methods place the 
source in the region between Seyfert and LINER galaxies ac- 
cording to emission line diagnostics dKauff^mann & HeckmanI 
2004). Due to large uncertainties in the [SII] lines fluxes, we can 
not use such measurements for a reliable input in the diagnostics 
for the source classification. 

In the X-ray regime, NGC 3599 was detected in a high state 
by XMM-Newton with respect to the RASS having a 0.2-2 keV 
EPIC-pn vs PSPC-2cr upper limit flux ratio of 88. In addition, 
the source also lies in the field of view of two ROSAT PSPC 
pointed observations and corresponding upper limits imply flux 
ratios of 297 and 100 respectively. Both XMM and ROSAT 0.2- 
2 keV fluxes were estimated from count rates using PIMMS as- 
suming a typical near-maximum tidal event spectrum, i.e. a black 
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Table 1. Spectral fit parameters (including 90% confident errors) 
of the EPIC-pn slew spectrum of NGC 3599. 
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0) 



Model 


Tx 


kT (keV) 


xHA.o.t.) 


wabs(bb) 


(...) 


09'i+o'"''* 


7.6(7) 


wabs(po) 


3 33+'-25 


(...) 


6.1(7) 


wabs(po+bb) 


2-86!ii? 


0.009!°:°°] 


4.2(5) 


wabs(po+MEKAL) 


^•^-^-0.64 


0.21!^% 


5.6(5) 



body model with kT—O.Ql keV modified by Galactic absorption 
(Komossa 2002). The unabsorbed X-ray luminosity in the soft 
0.2-2 keV energy band at the time of the slew observation was 
calculated to be 5.1(+0.7)xl0'*' erg s"' . As the source lies within 
the Virgocentric flow the luminosity distance is not a relevant 
parameter for the ca lculation of the X -ray luminosity, instead, a 
distance of 23 Mpc (iLauer et al.l l2007b) has been used. 

We have analysed the slew observation coiTesponding to the 
detection of NGC 3599. The object does not suffer from pile-up 
problems so the pn spectrum was extracted with patterns 0^. 
The source spectrum has been extracted from a circle of radius 
60 arcsec centred on the object. An adjacent region of the same 
radius has been used for the background. As software to produce 
response matrices for slew sources is not available, we have cre- 
ated the ancillary file using an ad hoc method of averaging the 
effective area along the track of the slew (see Read et al. 2008, 
for a full explanation of the technique). For the redistribution 
matrix we used the canned response file epnJf20jdY6jv6.9.rmf, 
appropriate for the average position of the source in the detec- 
tor. Due to the low signal-to-noise of the object, a binning of 
5 counts per spectral channel has been used for the spectral fit- 
ting. This is enough to perform a quantitative analysis although 
not statistically significant. Therefore, we can only attempt to de- 
rive a rough spectral shape of the object applying simple spectral 
models (Table [1]). 

As a result from the spectral analysis using XSPEC 11.3.2, 
a good fit to the data (x^^^^^-^^) is provided by a black body 
model of kT=95^l eV modified by Galactic absorption (Fig.|2l). 
This temperature is consistent (within the errors) with measure- 
ments from pr evious candidat es although at the upper range of 
temperatures dKomossal |2002|) . The corresponding 0.2-2 keV 
X-ray luminosity from the XSPEC fit is 5.5 x 10"^' erg s '. 
Using a power-law model, an equally good fit is achieved with 
rx=3.3+[ j2. A range of spectral models for the high-state spec- 
trum have been used in PIMMS to estimate the flux at the time 
of the slew observation. For black body models with temperature 
within 0.04-0. 1 keV and power-law models with X-ray slopes of 
2.9-3.5, file luminosities lie in the range 4.4 x 10'*'-6.5 x 10"*' 
erg s~'. The median luminosity from the different spectral mod- 
els is obtained with the black body model at 0.07 keV, i.e. the 
typical temperature for a tidal disruption event. The extrapola- 
tion of the slew count rate in the 0.2-12 keV to calculate an 
upper limit for the 2-12 keV luminosity gives 2 x 10"*" erg s"' 
for the power-law model with F^^ 
the black body model with kT^Q.Q95 keV). 



=3.3 (and 7 x 10^^ erg s"' for 



3.2. SDSSJ1 32341. 97+482701. 3 

XMMSLl J132342.3H-482701 was detected in 2003-12-01 
with a count rate of 2.08 cts s~* in the slew total energy 
band. Its counterpart, SDSS J132341.97-H482701.3 (named as 
SDSS J 1323 henceforward), can be classified as a non-active 
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Fig. 2. XMM-Newton EPIC-pn slew X-ray spectra of NGC 3599 
fitted by the black body model plus Galactic absorption. 
Data/model residuals are presented at the bottom. 



galaxy as its optical spectrum only shows absorption line fea- 
tures from the atmospheres of individual stars. Prior to its de- 
tection as tidal disruption candidate, no reference regarding this 
source is present in the literature. 

Optical data from a pre-outburst observation (performed 
seven months before the corresponding slew exposure) were ob- 
tained from the Sloan Digital Sky Surve>{3 (SDSS) and modelled 
using the EL-ICA method. No emission lines could be detected 
and the continuum can be best fitted by a single stellar popu- 
lation of 6.25 Gy and a metallicity of Z=0.4 Z© (Fig. [3]i. The 
galaxy is at a redshift of z=0.0875 and its stellar velocity dis- 
persion has a value of cr=80(+10) km s"' as derived from the 
optical spectrum. As no emission lines that could bias the source 
classification appear in the spectrum, no additional information 
can arise from applying the CONFIT approach to the stellar sub- 
traction. 

The X-ray slew detection of the object showed a flux varia- 
tion in the 0.2-2 keV energy band of a factor 83 between the 
slew detection and its RASS upper limit. The redshift of the 
source implies a luminosity distance of 399 Mpc, that allowed 
us to derive a soft unabsorbed X-ray luminosity of 4.8 x 10"*"' 
erg s"' from the slew observation using the black body model 
for tidal disruption events near outburst. The slew observation 
contained too few counts, all of them in the 0.2-2 keV energy 
band, to allow spectral fitting. 



4. Follow-up observations 

4.1. X-ray observations 

4.1.1. XM/W-A/ewfon observations 

XMM-Newton observations of 5 ks of both NGC 3599 and 
SDSS J1323 were performed roughly two years after the cor- 
responding slew exposures and then analysed with SAS 7.1.0 
(see Table |2] for details of the observations). All EPIC-pn and 
MOS observations were taken in Full Frame (FF) mode with 
the medium filter in place. High background periods produced 
by intense soft proton fluxes were rejected during the data re- 
duction process by defining good time intervals within the com- 
plete observation window. Source photons were extracted from 
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Fig. 3. Detailed fits to the optical spectrum of SDSS J 132341. 97+482701. 3. Flux calibrated data from the SDSS server are presented 
in black. The best modelling to the continuum as derived from the EL-ICA approach is shown in red and corresponds to a stellar 
population of 6.25 Gy and metallicity Z=0.4 Zq. The residual spectrum resulted from the subtraction of the best fit model to the 
observed spectrum is shown in blue. 



Table 2. Summary of the X-ray follow-up observations. 



Source name Mission/Instruinent ObsID 



Exposure 

(ks) 



XMM-Newkm /EPIC 2006-06-23 5 

Swift/XKT 2006-12-01 5 



NGC 3599 



XMM-Newkm /EPIC 2006-07-15 5 

Swift/XKT 2007-01-18 5 



SDSS J1323 



tion, hence all spectral fits were performed assuming Galactic 
foreground absorption values (A^h=1-42 x 10^° and 1.17 x 10^° 
cm 2 for NGC 35 99 and SDSS J 1323 respect ively) inferred from 
HI observations dDickev & Lockmanlll990l) . For both objects 
the hard X-ray behaviour could not be perfectly constrained by 
model fitting as a single data bin covering a wide energy range 
is present in that area of the spectra. Detailed analysis of the in- 
dividual objects will be described below and spectral parameters 
have been compiled in Table |3] 



circular regions of 30 arcsec radius centred on the object po- 
sitions. Circular source-free regions on the same chip and ra- 
dius of 60 arcsec have been used to characterise the background. 
The targets do not suffer from pile-up problems, so pn spectra 
were extracted with patterns to 4 and MOS spectra from 
to 12. Ancillary files and response matrices were generated by 
the arfgen and rmfgen SAS tasks respectively in order to con- 
vert the counts to physical units during the spectral analysis. All 
EPIC spectra were binned to oversample the instrumental reso- 
lution and to have no less than 20 counts in each background- 
subtracted spectral channel. Due to the low signal-to-noise ratio 
of the MOS spectra in both objects, the grouping resulted in at 
most three spectral points from each MOS observation. As a con- 
sequence, we only considered pn data for the source modelling 
performed using XSPEC. 

We initially performed spectral fits on both sources using a 
simple black body model with cold absorption to draw a com- 
parison with the high-state spectrum of a typical tidal disrup- 
tion event. In both sources the black-body model resulted in a 
poor fit. The post-flare spectra appear to have a roughly power- 
law shape. One-temperature thermal models often used to fit 
the spectra of normal galaxies (Raymond-Smith and MEKAL; 
dFab biano 1989)) also do not fit well, indicating that we continue 
to see the remnant of the outburst. From the spectral analysis, 
no improvement was achieved when including intrinsic absorp- 



X-ray emission from NGC 3599 has been detected in both 
soft 0.2-2 keV and hard 2-12 keV energy bands with corre- 
sponding count rates of 0.115(+0.007) and 0.010(+0.003) cts 
s ' . In order to give a first estimation of the spectral shape we fit- 
ted a power-law model and obtained a reasonable fit ix^^A-^-^^) 
with slope of rx=3.02;^Q2j but showing considerable residuals 
especially above 3 keV (Fig. |4] left panel). Models including 
ionized absorbers or partial covers {absori or zpcfabs in XSPEC) 
as well as those combinations of power-law and black body or 
thermal component (both using abundances fixed to solar value 
and treated as free parameter for the MEKAL model) do not 
provide a significant improvement. In all cases the F-test gives 
a confidence level of less than 85% with respect to the single 
power-law fit. In addition, the reliability of this models is not 
demonstrated as 90% confident errors on the parameters can not 
be constrained within sensible values. In any case the presence 
of a steep power-law is always required. Some spectral features 
remained unfitted as shown by the model residuals but with the 
present data no further constraints can be derived in this respect. 
This is just an indication that the real model is more complex 
than the ones here proposed as several emission mechanisms 
might contribute to the observed X-ray spectra. The unabsorbed 
luminosity at the time of the XMM-Newton observation derived 
from the simple power-law model is 1.8;^q4 x lO'"^ erg s"', i.e. 
a factor 27 lower than during the slew exposure. By extrapolat- 
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Table 3. Spectral fit parameters of the EPIC-pn spectra of the sources. 90% confident errors are quoted when;t'^^^ <2.0 
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Fig. 4. XMM-Newton EPIC-pn X-ray spectra of NGC 3599 (left panel) and SDSS J132341.97H-482701.3 (right panel) fitted by a 
po'wer-la'w model plus Galactic absorption. Data/model residuals are presented at the bottom of each spectrum. 



ing the best fit power-law model to high energies, we obtained a 
2-12 keV luminosity of 0.2+;}[ x lO^o erg s"'. 

For SDSS J 1323 the simple power-law model modified by 
cold absorption yields a reasonable fit in the soft X-rays but 
leaves an apparent high energy excess (Fig.|4] right panel). This 
model can be best fit using a steep power-law index rx=3.37 and 
the resulting x^^^^ is 2.49. Confidence errors on the slope could 
not be derived due to the high value of the;i'^g^. The same models 
as for NGC 3599 have been used in the analysis of SDSS J1323. 
Although the derived xl^.^ is closer to unity using some of the 
complex models, unusually steep slopes are required and fitting 
parameters can not be constrained. No significant improvement 
was achieved with respect to the simple power-law model as con- 
firmed by the F-test. Given that the spectrum seems to present an 
absorption feature around 0.5 keV, we tried to fit it with a power- 
law plus edge model. A simple edge at an energy of 0.49+qJ2 
keV gives a good fit ix^^^ - 0.56) albeit with a very small num- 
ber of degrees of freedom. This is consistent with absorption by 
cold Oxygen at the redshift of the source. Although the same 
feature appears when grouping the data with different spectral 
binnings the low signal-to-noise precludes further interpretation. 
Count rates in the soft and hard X-ray bands are 0.022(+0.002) 
and 0.002(+0.001) cts s^' respectively. The inferred 0.2-2 keV 



X-ray luminosity as derived from the simple power-law model 
is 1-1 ^05 X 10"*^ erg s"', implying a decrement of a factor 41 in 
comparison with the previous slew observation. 



4.1 .2. Swift observations 

Follow-up observations of 5 ks have been performed with the X - 
ray Telescope (XRT) on board the Swift satellite (Gehre ls et al.l 
2004) within the fill-in time program (Table|2]i. Data from pho- 
ton counting mode were analysed using XSELECT adopting 
the standard 0-12 grade filtering. The observing mode that has 
been used operates in the energy range 0.2-10 keV and retains 
full imaging and spectroscopic resolution. Target positions have 
been determined using the xrtcentroid tool. Source photons have 
been extracted from circular regions of 2 1 .2 arcsec based on the 
count rat es of the source s (as inferred from the recommenda- 
tions in Evan s et al]|200 7^ centred on the objects. For the back- 
ground determination, annular uncontaminated regions around 
the sources have been selected containing around a hundred 
counts in order to minimise the effect of background fluctua- 
tions. 

Both objects turned out to be extremely faint at the time of 
the Swift observations and, with 23(±5) and 10(±3) counts for 



p. Esquej et al.: Evolution of tidal disruption candidates 



1000 



500 




k««AlU 



o ^/V»/%i«nHwv^^ 



M^w»ljVl^J^vlM1fl^v^^ 



***AW*> 



J 

4000 



\ 

4500 



\ \ 

5000 5500 

rest wavelength (A) 



\ 

6000 



6500 



Fig. 5. Detailed fits to the post-outburst optical spectrum of NGC 3599. Flux calibrated data are presented in black and the stellar 
population is shown in red. The subtracted spectrum in blue illustrate the galactic nuclear region. 



NGC 3599 and SDSS J1323 respectively, we can not apply spec- 
tral fits to model the data. In order to estimate the source fluxes 
we applied a count rate to flux conversion through PIMMS using 
the same spectral model derived from the previous XMM pointed 
observations (power-law with rx=3.02 and 3.37 for the different 
objects according to section 14. Lit to our simulations. This is 
the only reasonable assumption that we can make because the 
behaviour of the post-outburst spectra is a priori unclear. With 
this hypothesis the source flux of NGC 3599 declined by a fur- 
ther factor of two in comparison with the previous X-ray ex- 
posure performed six months earlier. SDSS J 1323 is consistent 
with constant luminosity within the errors. 

4.2. Optical observations 

Long-slit optical observations of NGC 3599 and SDSS J 1323 
have been performed in order to obtain post-outburst spectral 
shapes and to search for possible signatures of the tidal disrup- 
tion events. 

Optical spectra have been obtained in December 2007 
with the ALFOSC spectrograph on the 2.5 m Nordic Optical 
Telescope (NOT) in La Palma using the grism 7, that covers a 
total range from 3850 to 6850 A with an instrumental disper- 
sion of 1.5 A per pixel. Exposure times were 1200 s and 1800 
s for NGC 3599 and SDSS J 1323 respectively. The calibration 
lamp and flat-field exposures were taken before or after any sci- 
ence exposure at the position of the target in order to keep con- 
trol of possible dispersion changes with the telescope position. 
Data reduction was performed using the IRApQ and STARLINK 
software packages. Initial reduction of the CCD frames involved 
bias substraction, flat-field correction and cosmic rays removal. 
In order to obtain the nuclear spectra, central apertures on the ob- 
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jects have been extracted. Wavelength calibration has been per- 
formed using ThAr lamp exposures. The wavelength calibration 
accuracy reached was always better than 0.5 and 1 A in the resid- 
uals for NGC 3599 and SDSS J 1323 respectively. 

From the flux calibrated nuclear spectrum of NGC 3599, we 
obtain [N II] 6584 A/Ha=l.l(+ 0.1). This emission Une ratio 
stays constant with respect to the pre-outburst observation. The 
[O III] 5007 A/Hy8 ratio difi'ers from that of the pre-flare expo- 
sure. This is due to calibration issues in the optical data and not 
from real variability of the lines as could be stated from a poste- 
rior observation of the object (see below). The current spectrum 
does not exhibit any dramatic emission-line response, such as 
prominent broad Balmer lines or strong high-ionization lines as 
it has been recently discovere d in the tidal dis ruption candidate 
SDSS J095209.56+214313.3 (lKomossaetal.1 12008). The post- 
outburst spectrum of SDSS J1323 showed no emission lines nor 
apparent variation since the spectrum performed before the dis- 
ruption event. Due to the faintness of the source and the rela- 
tively short exposure time of the observation, a reliable flux cal- 
ibration could not be achieved. However, this is not required for 
our aims as it has been confirmed that the source does not show 
at present any evident optical signature of the flaring event. 

A new optical observation of NGC 3599 was performed in 
April 2008 with the IDS spectrograph on the 2.5 m Isaac Newton 
Telescope (INT) in La Palma. Given the strong underlying host 
galaxy absorption features, uncertainties in existing spectra are 
too large to assess small-amplitude variability of the emission 
lines. This observation was aimed to obtain for the first time a 
high quality spectrum of the source. Three exposures of 1800 s 
of the object were performed using the R300V grating and then 
analysed following the same procedure outlined for the NOT 
observations. A S/N~90 was achieved in the final spectrum. 
The linear fit applied between pixel number and wavelength 
for the selected spectral lines gave a calibration error less than 
0.2 A. A set of 14 star exposures were observed and analysed 
in order to optimise the flux calibration process. Sensitivity 
functions of all standards were derived and combined into a 



p. Esquej et al.: Evolution of tidal disruption candidates 



unique solution using an iterative process that rejected those 
observations whose sensitivity behaviours deviated more than 
two percent from the mean calibration function. Gaussian 
profiles were fitted to the emission lines of the resulting 
nuclear spectrum of NGC 3599 after the starlight component 
subtraction (Fig. |5]l. The following ratios were obtained: 
[O III] 5007 A/Hy3=3.38(+ 0.91), [N II] 6584 A/Hq'=0.99(± 
0.07), [S II] 6717+6734 Aftia^0.3l(± 0.05), 
[O I] 6300 A/Hq'=0.11(+ 0.03). The source is located 
in the transition region between Seyfert and LINERS 
dKauffma nn & Heckman 2004) and, presenting a Lh„ = 10^^ 
erg s~ ', can be considered as a low luminosity AGN dHo et alj 
fT997h . 

5. X-ray light curves of the flare events 

The X-ray flare begins when the most bound material returns 
to pericentre after completing one orbit. Assuming that the 
mass distribution of the debris is near ly uniformly distributed as 
derived fro m numerical simulations ( JEvans & Kochaneklll989l: 
lAval et al]|200Q) . the mass rate after one post-disruption orbit is 



M: 



t - fdis 

lyr 



-5/3 



(2) 



where tdisr is the time at which the disruption process occurred. 
Luminosities in the soft 0.2-2.0 keV energy band have been 
derived from all three X-ray source detections. Assuming that 
the debris evolution of our sources is represented by the fallback 
model of tidal disruption, data points in the light curve have been 
fitted with a (f - fdisr)"^^'' dechne law. These fits give 



Lx = 6.7(±1.2)xl0' 
and 



,40 



f-2003.59(+0.04)yr 



1-5/3 



1 yr 



ergs 



(3) 



event has been estimated for both objects. The associated total 
masses involved in the process have been derived by applying 



AM = 



Afi A£x 



(6) 



where e is the efficiency of converting mass to radiated energy 
(we used a value of e=0.1). It has been assumed that the bolo- 
metric correction is not very significant, which is a fair approx- 
imation if the spectrum is well represented by a black body as 
we expect near outburst. Assuming a black body model at 70 eV 
at the corresponding redshift of our sources, the bolometric cor- 
rection would imply less than a factor of two. These results are 
presented in Table |4l together with the black body radius of the 
emitting region (as derived in lLi et a02002l) expressed by 



R^ 






1/2 



(7) 



where /(, = T\,\,/T^ff and has been set to /c=3 for our calcula- 
tions. Although the luminosity behaviour is not well constrained 
in the vicinity of the flare onset, we can also perform our calcu- 
lations assuming tha t the beginnin g of the flare occurs at a time 
after the Af a; 0.021 (iLiet all 12002. using a 10*' Mq black hole 
and a solar-type star), after the disruption. In that case, the esti- 
mated parameters in Table |4] would increase by about an order 
of magnitude. 

Table 4. Estimated properties of the disruption events. Column 
1: source name. Column 2: energy released during the observed 
outburst. Column 3: accreted mass associated with the total 
event. Column 4: black body radius of the emitting region. 



Source 


A£x (erg) 


AM (Me) 


Rx (cm) 


NGC 3599 


7.1(±0.8)x lO'"* 


4.o(±o.4) X ur^ 


7.3(±2.7) X 10" 


SDSSJ1323 


7.6(±0.3) X lO'O 


4.2(±0.2) X 10-3 


6.8(±3.4) X 10'2 



8.1(+2.9)x 10"* 



f-2003.56(+0.10)yr 



1-5/3 



1 yr 



ergs 



(4) 



for NGC 3599 and SDSS J1323 respectively. Fig. |6] shows the 
evolution of the X-ray luminosity for both sources and corre- 
sponding fitted light curves according to the previous equations. 
Although for the slew detections the representative black body 
model ikT-Q.Qil keV) for tidal disruption events has been used 
for the luminosity calculation, error bars for these data points 
comprise the values obtained when allowing a range of mod- 
els for the outburst spectra (black body with temperature within 
0.04-0. 1 keV and power-law with X-ray slopes of 2.9-3.5). 

In order to estimate the total X-ray energy released during 
the event, the previous luminosities have to be integrated over 
the outburst duration via 



^E^ 



-r 



Lx{t)dt. 



(5) 



The flare, that begins when the most bound material returns to 
the central black hole after completing one orbit, occurred be- 
tween the disruption time previously calculated (see eqs. |3]and 
m and the slew detection. We will use the time of the slew obser- 
vations as conservative assumption for the beginning of the flare 
(lower limit of the integral in eq.|5]l in the following calculations. 
With this hypothesis, the total energy released during the entire 



Several scenarios can be argued in order to explain the low 
amount of accreted matter. First, it may imply that we have only 
seen a partial disruption event. Conversely, it can also be in- 
dicative of a strong encounter when the star deeply penetrates 
within the tidal radius leading to a tidal detonation event. In 
that latter case, the star would experience a strong pre-disruption 
compression triggering thermonuclea r reactions and being fol - 
lowed by the explosion of the star (e.g. lCarter & Luminetll982h . 
Therefore, only a small fraction of the debris would be ac- 
tually accreted and the r est of the material being expelled 
dBrassart & Luminetl2008h . The tidal radius (for a Sun-like star) 
as given by eq. O) and the Schwarzschild radius are 



Rj -7x10' 



i 



M, 



BH 



\ IO^Mg 



1/3 



M, 



\10''Mc 



cm, rs ~ 3x 10 I ^^^^^ | cm. 

(8) 

The possibility of a stong encounter would be reinforced, espe- 
cially for NGC 3599, as the derived radius of the emitting region 
lies inside the tidal radius for a 10^ Mq black hole. Nevertheless, 
further interpretation is out of scope as these values are simple 
order of magnitude estimates. 

It is demonstrated that X-ray and optical luminosi ties of 
early-type galaxies present a tight correlation (O' Sulliva n et al.l 
2001). Studies on NGC 5905, a previous case of tidal disrup- 
tion, showed that the optical luminosity of the source did not 
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Fig. 6. X-ray light curve of NGC 3599 (left panel) and SDSS J132341.97+482701.3 (right panel). In each plot, data points (90% 
confident errors are included) correspond to XMM-slew, XMM-pointed and Swift observations from left to right. Solid lines are 
the fitted curves predicted for the fallback stage and follow eqs. ([3]) and (|4]l respectively for both sources. Dashed lines show the 
disruption time (fdisr) as returned by the fitting. 



experience any significa nt variation due to the flaring event 
dKomossa & Badelll999l) . We will assume that this is also ap- 
plicable for our sources and we will use pre-outburst photomet- 
ric data for our calculations (no relevant low-state information 
is available for comparison). B band magnitudes were obtained 
from the NED database (for SDSS J 1 323 conversion factors have 
been applied to the Gunn magnitudes to derive the magnitude in 
the B Johnson filter). At the time of the Swift observation we 
derive ratios (log Lx/Lb)=30 and 32 (with Lg in units of solar 
luminosity) for NGC 3599 and SDSS J 1323 respectively. These 
values are beyond the observed limiting value for early-type 
galaxies and, therefore, we can argue that the X-ray emission 
of the sources is most likely related to the flare event. 

An alternative estimation to characterise the nature of the 
low state emission of our sources can be reported by using the 
correlation of [OIII] and hard X-ray luminosities (L2-10 kev vs 
L[om]) of Seyferts and LLAGN. L[oin] luminosities have been 
used to predict the expected har d X-ray luminositie s following 
eq. (1) in the recent estimation of iNetzer et al.l ( 12006 ) from their 
X-ray selected sample of sources. These values have been then 
compared with the X-ray luminosities in the 2-10 keV energy 
band obtained by extrapolating the best fit power-law model in 
the XMM-Newton pointed observations. As can be seen from 
Fig. m sources are not very well modelled in this energy range 
and data points lie above the model, implying that our estima- 
tion falls below the real value of the luminosity. No significant 
change in spectral shape due to the flaring events is present in 
the optical data of both sources, indicating that the we could use 
either the pre- or post-outburst value of the [OIII] emission. For 
our calculation we used the [OIII] value from the highest res- 
olution available spectra for both sources, i.e. the most recent 
optical observation of NGC 3599 and the pre-outburst observa- 
tion for SDSS J1323. For NGC 3599, the estimated 2-10 keV X- 
ray luminosity through the [OIII] correlation is 6 x 10^^ erg s '. 
This value is consistent (within the errors) with the 1.5 x lO-'^ 
erg s"' luminosity returned by the spectral fit from section [4.1. 11 
Although we can not constrain with this argument that the low- 
state X-ray emission is not AGN related, it is worth mention- 



ing that the equation used for predicting the X-ray luminosity 
was derived for luminosities down to L2-iokeV=10'*' erg s"'. As 
NGC 3599 presents a luminosity two orders of magnitude lower, 
this can conflict with the validity of the formula for our calcula- 
tion. The optical spectrum of SDSS 1323 does not present emis- 
sion line features. A single gaussian fit has been performed to 
the position of the [OIII] emission line and a 3cr upper limit on 
the L[oui] has been derived, with redshift and linewidth fixed to 
that determined from stellar absorption lines. In this case, the 
-^2-10 kev=3.4 X 10 ^^ erg s"' obtained from the relationship in 
iNetzer et all (l2006h is below the 4.3 x 10"*° erg s ' obtained from 
the spectral fit. As the low-state X-ray luminosity for SDSS 1323 
is higher than that predicted from the [OIII] emission, we have 
an extra argument to confirm that the low-state emission is flare 
related, rather that coming from a permanent AGN. 

6. Black hole mass estimates 

Due to evidence for the existence of a massive black hole lo- 
cated at the centres of most galaxies, a number of empirical re- 
lationships have been derived relating the black hole mass and 
different characteristics of the corresponding host galaxies. 

Among them is the Mbh - cr relationship, that is the correla- 
tion between the black hole mass and the stellar velocity disper - 
sion. Using the recent estimation from lFerrarese & FordI (l2005l) 

/ a- x4.86(±0.43) 

Mbh = 1.66(+0.24)x lO^Mo — -] (9) 

\ 200 km s ' / 

we derive a black hole mass for NGC 3599 of 1.3(+0.6)xlO^Mo 
using 0--13.2 km s"' (several velocity dispersion measurements 
exist for this galaxy and the median value provided by Hyperleda 
as been used in our calculation). For SDSS J 1323 the mass is 
2.2(+O.9)xlO''M0, with cr=80 km s"' as obtained from the anal- 
ysis of the optical spectrum (see section |3T2] i (no other measure- 
ment is available for comparison). 

Con versely, using the Mbh - cr estimate from iLauer et al] 
(l2007bh 
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Mbh = 1.95(+1.17)x 10** Mo 



200 km s" 



4.13(±0.32) 



(10) 



for 



we obtain 3.1(+l.O)xlO''M0 and 5.O(+1.4)xlO''M0 
NGC 3599 and SDSS J 1323 respectively. 

Other estimation arises from the relationship between the 
black hole mass and the V absolute magnit ude of the host galaxy . 
Using the correlation derived in eq. (6) of iLauer et al.l (l2007ah . 
black hole masses of 3.5(+O.9)xlO^M0 and 3.3(+0.9)xlO^Mo 
are obtained for NGC 3599 and SDSS J 1323 respectively. 

Black hole mass estimates from both Mbh - o" and Mbh - Ly 
approaches differ by one order of magnitude approximately, due 
to uncertainties in the measurements and scatter in the correla- 
tions. Nevertheless, results from both calculations lie within the 
mass range expected for objects experiencing a tidal disruption. 

Using the black hole masses derived from the different ap- 
proaches, the observed high-state luminosity results in a frac- 
tion of the Eddington luminosity. That factor is up to 2% for 
NGC 3599 and 17% for SDSS J1323 (using the smallest black 
hole mass estimation). Therefore, at the time of the slew observa- 
tion the emission was sub-Eddington, especially for NGC 3599. 

7. Tidal disruption rate from the slew survey 

The identification of tidal disruption events discussed here was 
enabled by the existence of two large area sensitive X-ray sur- 
veys performed at different epochs. We can only attempt to cal- 
culate a first order estimate of the outburst rate as a number of 
assumptions have to be taken into account. 

In order to determine the characteristic outburst flux limit 
to which our survey is complete, we first derived the nominal 
RASS flux by computing 2cr upper limit count rates of 5000 
sources randomly distributed over the whole sky, as we assumed 
that the outbursting galaxies are isotropically distributed. These 
rates in the 0.1-2.4 keV energy band were converted to 0.2-2 
keV unabsorbed fluxes through PIMMS assuming the typical 
spectral model expected near the peak of the tidal event, i.e. a 
black body model with temperature 0.07 ke V modified by galac- 
tic absorptiorQ. The flux distribution is widely spread due to 
the different exposure times of the RASS observations and the 
different absorption column densities of the sources. In order 
to avoid high galactic absorption and confusing sources in the 
galactic plane, in the remaining calculation we will not consider 
the region with Galactic latitude |b| < 20 deg. From the flux his- 
togram, we derived a median flux of 2.5 x lO"'-' erg s"' cm"^. 
We will adopt this number as the characteristic unabsorbed 
flux limit for the RASS, which is in fair agreement with the 
3x 10"'^ erg s"' cm"^ typically quoted jBrandt & Hasingen 
l2005h . Therefore, if we consider as outbursting events those 
sources which varied by at least a factor of 20 with respect to 
the RASS, the characteristic flux to which our survey is com- 
plete is 5.08 X 10"'^ erg s"' cm~^. The on-axis flux limit in the 
soft band for slew observations applying the black body model is 
4.33 X lO"'-' erg s"' cm"^ and raises to 4.33 x 10"'^ erg s"' cm"^ 
for the minimum exposure times, hence we reach the required 
sensitivity for completion derived from the RASS picture in the 
whole slew coverage. Tidal disruption events near maximum are 
expected to emit a luminosity up to 10"*^ - 10'*^ erg s"' de- 
pending on the black hole mass. Using an outburst luminosity 




0.00 



0.02 



0.04 0.06 

redshift 



0.08 



Fig. 7. Observable flare duration depending on redshift. 



of 10"^'* erg s"' (the median value for our two candidates) and 
the flux limit above calculated we derived that the completeness 
distance of our survey is 406 Mpc (z=0.089) that, in terms of 
volume, corresponds to 2.8 x 10^ Mpc^. 

In order to derive the tidal disruption rate we will have to 
consider the time interval in which the flare, that starts when the 
most bound material returns to pericentre, is visible before the 
source has declined to non-detectability. To estimate this figure, 
we will assume that for nearby sources the galaxy luminosity 
peaks at 10"*"^ erg s"', follows the r^l^ law and then declines to 
its quiescent value in three years after outburst. Whatever the 
morphologycal type, the lumin osity of a n ormal galaxy lies in 
the range 10^*^-10^^ erg s"' (Fab bianol 19891) . so we presupposed 
a dormant luminosity of 10"**^ erg s"'. The next step will be to 
simulate that curve including the redshift dependence by means 
of the flux decrease due to the geometric dilution factor, as the 
unabsorbed X-ray luminosity depends on the squared luminosity 
distance (Lx = Fx^-nD^). With these hypotheses, all of them 
in fair agreement with theory, we derive that for sources up to 
redshift 0.0025 we will see the events at least for one year before 



declining below our limiting flux of 5.08 x 10" erg s" cm" . 
The observable duration of the flare depending on redshift and 
estimated under these premises is presented in Fig.|2] 

To obtain the number of events we detect per year per unit 
volume we will divide the number of tidal candidates by 



f 

Jo 



A(r)t(r)dr 



(11) 



^ Galactic absorption for each source position was calculated using 
the Nii tool provided by PIMMS. 



where ^max=406 Mpc as derived before, A(r) represents the area 
covered (including the 16.5% fraction of the sky surveyed by 
slew obsevations in the XMMSLl_DeltalQ with |b| < 20 deg), 
and t(r) is the observable flare duration as observed in Fig.|2] As 
we detected two tidal diruption candidates, the outburst rate we 
obtain is 5.4x10"'' yr"' Mpc"^ 

In order to estimate the rate of events per galaxy we need to 
take into account the total galaxy space density. Extensive work 
has been done in deriving luminosity functions of galaxies de- 
pending on their morphology, as different galaxy types require a 
distinct modelling using either Gaussian or Schechter functions. 
That topic is beyond the scope of the present paper and we will 
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simply use the derived parameters from lde LapparentI (l2003h for 
our calculation. In the same way as in other attempts of calculat- 
ing the tidal disruption rate, i.e. including all E, SO and elliptical 
galaxies in the B band for the space density calculation, we de- 
rived a total space density of 2.3 x 10"^ Mpc"^. Therefore, the 
resulting tidal disruption rate is 2.3 x lO""* galaxy"' yr"'. 

In order to validate our result, we have compared our tidal 
disruption rate with that derived from theoretical pred i ctions . 
The tidal disruption rate as estimated in I Wang & Merrill (|2004|) 
can be expressed by 



2.4- 



fc 



1.6- 



1.4- 



A^ 



7.1xlO"V"' 



7/2 



Mr 



-1/3 



-1/4 



R*^ 

70kms-i/ llO^Moj [MsJ Uo 

(12) 

cr being the velocity dispersion of the host galaxy, and m* and 
Ri, the corresponding mass and radius of the disrupted star. This 
equation was normalised to reproduce the tidal disruption rate 
for cr-7Q km s ' and a 10* Mq black hole. The dependence 
on the stellar mass and radius is weak as the tidal disruption 
rate is expected to be dominated by subsolar stars (which ful- 
fill the relation R*/Rq ~ mi,/MQ). We will assume solar mass 
and radius for simplicity. Merging eqs. (|9]l and ( fTSl l. the unique 
dependence of the tidal rate on the black hole mass becomes ex- 
plicit. Considering a 10* Mq black hole, the consumption rate is 
T.OxlO"'* galaxy yr"'. Note that eq. ( fT2] ) only provides a good 
parameterisation for power-law galaxies with F > 0.9, overesti- 
mating the tidal disruption rate by up to an order of magnitude 
for galaxies with shallower power-law and core brightness pro- 
files. Although a number of simplifications have been included 
in our estimations, our predicted tidal disruption rate as derived 
from the observations agrees with that from the theoretical ap- 
proach. 

For completeness, although not very important as we are 
dealing with nearby objects, we have estimated the dimming ef- 
fect on the flux determination due to the redshifted energy band, 
i.e. the A'-correction, that in our case can be expressed by 



A'n?- 



0.2-2.0 keV - 



lo.'ify'^'' 



/o 



2.0(l+z) 
'0.2(1 +<•;) J''^'^ 



(13) 



where fy is the unabsorbed spectral flux density. X-ray spec- 
tra have been simulated in XSPEC to derive that bandpass 
correction for a redshifted black body model (zblackbody) at 
kT-lQ eV using both EPIC-pn and ROSAT response files. The 
galactic absorption has been taken as A^h=4.03 x 10^" cm"^, 
which is the median value for the random positions in the sky 
calculated at the beginning of this section. This effect, that is sig- 
nificantly smaller than the geometric dilution factor, produces a 
decrease up to 10% for the XMM and 27% for the ROSAT fluxes 
at our redshift limit (Fig.[8]i. 

We have performed our calculations presupposing that in- 
trinsic absorption does not prevent the detection of the flaring 
events. As discused in Sembav & West (1993), this eff'ect can 
by included by estimating the rate of non edge-on spiral galax- 
ies in our sample. Due to the large uncertainty existent in the 
derivation of that number, the former calculation did not include 
this argument. As a rough approximation, assuming a fraction 
of 60% for object s in which intrinsic absorption does not play 
an important role (iDonlev et al.ll2002l) . the tidal disruption rate 
would increase up to 3.9 x 10"^ yr"' Mpc""'. 




Fig. 8. K-coiTection. Solid line: XMM-Newton EPIC-pn. Dashed 
line: ROSAT PSPC. Dotted line: limiting redshift for complete- 
ness in our calculation. 

8. Summary and conclusions 

We have presented X-ray and optical follow-up observations of 
two tidal disruption candidates detected by XMM-Newton dur- 
ing slew observations. 

The softness of our sources near outburst confirms that tidal 
disruption events are certainly very soft. Both NGC 3599 and 
SDSS J1323 were not detected in the hard 2-12 KeV slew en- 
ergy band close to the peak of the outburst. As all previous X- 
ray candidates were detected with ROSAT, the spectral softness 
near outburst could have been biased due to selection effects. 
However, this is unlikely because they were very soft even within 
the ROSAT 0. 1-2.4 keV spectral coverage. 

Although tidal disruption phenomenon is the most plausible 
scenario in order to describe the events here discussed, extreme 
variability of a possible ultraluminous X-ray source (ULX) lo- 
cated within NGC 3599 can not be safely excluded with the 
present data due to the relatively low luminosity of the source 
near outburst. Nevertheless, a ULX interpretation appears to be 
a priori improbable as this class of objects are characterised by 
X-ray luminosities in the range 10-'''-10'*'^ erg s"', flux variability 
factors of 2-3 and a power-l aw spectral shape wit h X-ray photon 
index between 1.6 and 1.8 tColbert & Ptakll2002l) . 

XMM-Newton and Swift observations of the two tidal dis- 
ruption candidates described in this paper have been performed 
in order to monitor the evolution of the objects. These are the 
nearest hard (> 2 keV) X-ray follow-up observations to an X- 
ray high state detection ever performed. We can now assess the 
existence of hard X-ray photons within three years after the ob- 
served high state, although a spectral hardening can not be quan- 
tified as this is highly dependent of the spectrum assumed near 
outburst. Prior to this work, a single reliable hard band mod- 
elling of the phenomena was only achieved ten years after the 
maximum (Komossaet al. 2004). No absorption excess to the 
Galactic value is required in order to fit the spectra from the 
present data, although there is some evidence for an edge in the 
SDSS J1323 spectrum. With power-law slopes of Fx=3 and 3.5, 
the spectra are still very soft compared with normal AGN. We 
tentatively interpret the hard emission as coming from an accre- 
tion disk although in the case of NGC 3599 emission from the 
galaxy itself cannot be excluded. 
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The X-ray flux is seen to decrease between the XMM- 
Newton-slew and Swift follow-up observations by factors of 30 
and 40 respectively for NGC 3599 and SDSS J1323. The light 
curves plausibly follow the expected r^^^ curve. Post-flare opti- 
cal spectra of both sources did not reveal any dramatic emission- 
line response due to the flaring event. 

After detailed optical spectroscopic analysis, NGC 3599 
turned out to be misclassified as a inactive galaxy and appears 
to harbor a dwarf active nucleus. This fact is not unexpected 
considering that the maj ority of the nearb y galaxies host weakly 
active AGN-like cores JHo et alj 1200 lb . This interpretatation 
was also inferred for the previous tidal disruption candidate 
NGC 5905 due to the detection of faint HST-hnsed post-flare 
high ionization lines (iGezari et al.l l2003h . Nevertheless, given 
that there is no HST spectrum of NGC 5905 prior to the out- 
burst, the emission lines detected in this source could also have 
been excited by the flare itself as seen in a very recent candidate 
tKomossa et al...2008.) . 

The tidal disruption rate as derived from slew observations is 
2.3 X 10""* galaxy"' yr"', which agrees with previous theoretical 
estimations, and further endorses the view that these X-ray flares 
are the result of the tidal disruption by a SMBH. 

The detection of these phenomena represents a unique 
tool to unveil the dormant black hole population at the cen- 
tres of the flaring galaxies. A complete census of the relic 
SMBHs is fundamental to constrain the accretion history of 
the Universe. It has been suggested that tidal disruption events 
could fuel the majority of low luminosity AGN (LLAGN) and 
perhaps significantly con tribute to their luminosity function 
JMilosavlievic et all |2006|) . More flare events are expected to 
emerge in the near future based on the data collected during 
XMM-Newton slew observations besides future X-ray surveys 
Uke that planned with eROSlTA. 
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